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Abstract. The first order nonlinear ODE ip(t) +smip(t) = q(t),q(t) = B + Acosujt, 
where A, B,u> are real constants, is considered, the transformation converting it to 
a second order linear homogeneous ODE with polynoimial coefficients is found. The 
latter is identified as a particular case of the double confluent Heun equation. The 
series of algebraic constraints on the constant parameters is found whose fulfillment 
leads to the existance of solutions representable through polynomials in explicit form. 
These polynomials are found to constitute the orthogonal normalizable system. 



Nowadays, electronic devices based on the Josephson effect in superconductors and, 
in particular, Josephson junctions (JJ) pQ play the important role in the measurement 
technique, serving, in particular, the core element of the modern voltage standards 
The application needs lead to the growing importance of the theoretical and 
mathematical modelling of JJ properties. One of the theoretical tools commonly used 
for this purpose is the RS J (Resistively Shunted Junction) model Ej which applies, 
in the case of overdamped J J [S], the ODE 



Here the (real valued) function <p(t) called the phase and describing J J state is unknown 
while q(t) representing the external impact to J J (the appropriately normalized bias 
current supplied by an external source) is assumed to be given. The dot denotes the 
derivative with respect to the variable t (the appropriately normalized current time). 

The goal of the present notes is the discussion of some results concerning the 
equation and its solutions in the particular case of harmonic q(t) most important 
for applications. Thus we assume, without loss of generality, 



where A, B,u are some constants subject to the condition A ^ ^ w in order to 
eliminate trivial situations. 

As it was first noted by V.M. Buchstaber (see [H]), the first order nonlinear ODE 
is equivalent, for arbitrary q(t), to the system of two linear ODEs 



ip(t) + simpit) = q(t). 



(1) 



q(t) = B + A cos cut, 



(2) 



2x(t) 

m) 



x{t) + q(t)y(t) , 
[q{t)x{t)+y{t)\, 



(3) 
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where x(t),y(t) are the new unknowns. They are connected with tp(t) by the equation 

Thus ip is twice the phase of the complex quantity x — iy. 

In the case of q(t) defined by Eq. (|5|L the replacing of the independent real variable 
t with the complex valued variable 

z = exp(iu;t), (5) 
translates Eqs. (JHJ) to the equations 

4iuz x(z)' = 2x{z) + [2B + A (z + z' 1 )] y(z), 
Aiuz y{z)' = - [2B + A (z + z' 1 )] x{z) - 2y{z). 
Here the separate notations for the unknowns x, y considered as functions of the complex 
variable z, x(z) = x(t),y(z) = y(t), are employed, the prime denotes the derivative with 
respect to z. Multiplying them by z, the equations with polynomial coefficients result. 

One may extend the original meaning of x, y, treating them as analytic functions 
of the complex variable z which satisfy Eq. © everywhere in the complex plane except 
singular points. Then, obviously, a non-zero real or imaginary parts of any pair of 
such functions which prove smooth on some segment of the unit circle in C yield a real 
solution of Eqs. (|HJ) and, consequently, a (real) solution of Eqs. ([H),©- 

Now let us consider the following transformation replacing the unknowns x, y with 
the functions v = v(z),v = v(z) in accordance with equations 

v= i z~^ exp ( — (— z + z^ 1 ) ) (x — iy), 

(A < 7) 

v = (2ujz) 1 z 2 " exp I — (— z + z x ) j (x + iy). 

It is easy to show that the fulfilment of Eqs. (jBJ) is equivalent to the condition of the 
vanishing of the following two expressions 



v and z 2 v + —(z 2 + !)+(— + 1) z)v + —rv (8) 



\2uj y \u J J 4cu 2 

that, in turn, is equivalent to the fulfilment of the second order linear ODE with 
polynomial coefficients 



0. (9) 



d 2 z \2lo \u J J dz Auj 2 

The latter has the only two singular points, z = and z = oo. 
The Mobius transformations 

C = ^, (10) 
z — a 

where ( is the new independent complex variable and a is an arbitrary non-zero complex 
number, leads to representations of (jOJ) more clearly revealing the symmetry in the roles 
of its singular points. In particular, in the case a = i, one gets the equation 
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which proves to be is a particular instance of the double confluent Heun equation 
(DCHE) as it is given in [7], Eq. (4.5.11). It arises for the following set of parameters 
employed in a = 0, c = — (Bu -1 + l),t = iA(2u)~ 1 , A = (2iuM) _1 . It is also worth 
noting that the canonical DCHE representation (Eq. (4.5.1) in jjj) 



. , v -2 2 + c2; + t)^^ + (-a2; + A)?/(2;) = 

dz z dz 

results from © after the argument rescaling z — > 2iuA~ l z and corresponds to the 
parameters a = 0, c = 5c<j _1 + 1, £ = — (\Auj~ 1 ) 2 , A = ^o~ 2 . 

At the same time the most elegant form of the Mobius-transformed Eq. Q results 
with a — 1 in which case one gets 



d( dC & J d C 



0(12) 



(see [8 ). Here the singular points ( = ±1 are just the images of the only singular 
points z = 0, oo of Eq. (J0J). Hence there are no more singular points for Eq. (|12j) and, in 
particular, the point ( = oo is regular. It is also easy to show that the only effect caused 
by the transition to the reciprocal variable, ( — > is the inversing of the sign of the 
parameter A in (|12p. Another £ transformation which obviously preserves the form of 
Eq. (|T2"|) is the reflection ( — > — £. It leads to the combined inversion of parameter signs 
A — > —A, B —> —B. Thus, the signs of the both parameters A and B are irrelevant 
in the sense they can be (independently) reversed by means of transformations of the 
independent variable £. 

Yet another useful representation of Eq. © results from the the replacing of the 
unknown v by the new unknown P = P(z) by means of the substitution 

v = exp ( — — z J P. (13) 



2uj 

It is convenient to represent the resulting equation as follows 

z(zP' - nP)' - \iz{zP' - nP) + (fi - z)P' + \P = 0, (14) 

where the following new constant parameters n, fi, A replacing the equivalent triplet 
A, B, u in accordance with definitions 

f B A A I - A 2 1 2 

n= ~U J'^aJ' A = ^ r= (2^ 5 ~ A4 (15) 

are utilized. 

The representation (|14j) enables one to conjecture that if the parameter n assumes a 
non-negative integer value then this equation may admit polynomial solutions. Indeed, 
for n = 0, a constant is its solution for arbitrary u>, provided the constraint A = is 
additionally obeyed. Further, for a positive integer n, let us apply the ansatz 

n 

p=p n = j2 akzk > ( 16 ) 

fc=0 
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where the constant coefficients a k (dependent also on n) have to be determined. 
Substituting P n into (JHJ), one gets the following system of n + 1 linear homogeneous 
equations 

= Aa + fJLdi, (17) 
= n{n — k + l)ajt_i + (A — k(n — k + l))a k + fi(k + l)afc+i for k = 1, . . . , n — 1, (18) 
= /ia n _i + (A - n)a n , (19) 

for n = 1 the subsystem ()18|) being void. It admits a nontrivial solution if and only 
if the determinant A n (A, fi) of the following (n + 1) x (n + 1) dimensional 3-diagonal 
matrix 



(20) 



/ A n- 1 

/it • n A — 1 • n fi ■ 2 

n ■ (n - 1) A — 2 • (n — 1) ^ • 3 

|ii • (n - 2) A — 3 • (n — 2) 







V o 



A - (n - 3) • 4 ii ■ (n - 2) 

• 3 A - (n - 2) • 3 fj, ■ (n - 1) 

fi -2 A - (n - 1) • 2 n ■ n 

X-n-1 I 



vanishes. 



A n (A, //) = det$ = 



(21) 



is the algebraic equation, of degree n + 1 in A, constraining the parameters /i, A. Its 
fulfilment is the necessary and sufficient condition for the existing of polynomial solutions 
of Eq. (JHJ). 

If (J2IJ) is fulfilled, the polynomial coefficients a k can be calculated as follows. Having 
introduced their appropriately scaled ratios R k by means of the definition 



^ ak ~ l ^ k = a n ^ k - n ^±^- J] i? J forfc = l,2,...,n-l,(22) 

j=k+i 



k a k 



T(k + 



Eqs. (JTHJ) are equivalent to the recurrence relations 



A 



+ 



k = 1,2,... 



kik — n — 1) ' k{k — n — l)R k+ i ' ' ' ^ 

specifying, generally speaking, a continued fraction. However, for positive integer n, it 
is truncated at the (n — l)'th step in view of Eq. fi"9|) which is equivalent to the equation 

R n = l- \/n. (24) 

Then, making use of Eqs. (J2H1) , (J2H) , one may determine, step by step, all R k . In turn, 
they determine coefficients a k up to an arbitrary common factor. 

Further, let us note that Eq. (j22j) can be represented in the following matrix form 



Rk 
1 



{Z k R k 



Mi 



Rk+l 
1 



(25) 
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where Z k 




(26) 



Eq. ()25|) can also be interpreted as a linear map on the projective vector space of 2- 
element columns defined up to a nonzero factor. Iterating it and making use of ()24jl . 

one gets 



Rk 




1 





'n-1 



n ^ 



n M , 



X 



j — k 



n 

1 



(27) 



Here the symbol ^ denotes the product of matrices, where the factors corresponding to 
larger indices j are situated at right with respect to the lower index ones. 

Eqs. (J23Jl - (J27)l determine the set of coefficients a k and, then, the polynomial P n (z) 
for arbitrary A,/i irrespectively of the fulfillment of Eq. (|2ip. However, if the latter is 
not satisfied, such a polynomial cannot be a solution of Eq. Q14j) . Then some of Eqs. 
(fT7]) - (|T9")) have to be not satisfied as well. Since (fTHj) . (fT^ j) are automatically obeyed by 
the very meaning of Eqs. (|2*3|) - (f2*Tj) . it is the fulfillment of Eq. (jl4j) which is equivalent 
to the equation 

Ri = - £ (28) 



which is, in turn, connected with fulfilment of Eq. (|17|). Accordingly, if i?i is considered 
as the result of calculation with the help of the formula (J27)) for k = 1, the equation 



HI 

A 



n-1 



X 



1-^ 







(29) 



arises as the necessary condition of the fulfilment of Eqs. (|TTj) -()19J ) . In view of the close 
resemblance of algebraic structures of Eqs. (|2Tj) and (J2T?|) . it is then natural to suppose 
that the left-hand-side expression of Eq. (|2*9*j) is intimately connected with A n (A,/i). It 
is, indeed, the case, and the following representation of A n = det $ through the product 
of a finite set of 2 x 2 matrices M.; (|2*B]) takes place: 



A n (/x, A) = - 



A 
/' 2 



-j 

T n-1 

X 



Il M i 



3=1 



n 
n 



A 



(30) 



Similarly, there is the following explicit representation of the polynomial coefficients a k 
through the analogous matrix products: 



a k 



k(n + !-&;)! 



T n-1 

X 



X 



j—k 



n — A 



1. 



1. (31) 



Although, formally, it does not cover the case k — 0, a minor modification allows to 
compute do as well (one has to replace in the multipliers the integer parameter k with 
k + e, where e is a small real number, to carry out computation with k = 0, and to pass 
to the limit e in the result). 
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The polynomials P n reveal a remarkable symmetry concerning the coefficients in 
front of the "small" and "large" z powers. It can be discovered considering the following 
polynomial constructed from P n : 

P n (z) = z^z- 1 ) - tiPniz- 1 )). (32) 

A straightforward computation shows that it obeys Eq. (jl4j) if and only if P n (z) does. 
However, the polynomial solution of Eq. (fTlj) is unique up to a normalization. Indeed, 
the second solution linearly independent with P n (z) admits the following representation 
in quadratures: 

Qn = Pn J z n exp (/i (z + z- 1 )) P~ 2 d z (33) 

(the associated functions) . As opposed to P n , it is obviously singular in the point z = 0. 
Thus, having expanded P n through the basis P n ,Q n , it cannot involve any fraction of Q n 
in its 'content' and, thus, z n (P' n (z~ x ) — [iP n (z~ 1 )) oc P n (z). The constant proportionality 
coefficient can be fixed evaluating the equations above at the point z = 1. In this way 
one may obtain 

P' n {z) - fiP n (z) = e(2u)- 1 z n P n (z- 1 ), (34) 

where e 2 — 1. 

It is therefore shown that the fulfilment of Eq. (|3*4*j) is the necessary condition for 
the polynomial P n (z) of degree n to satisfy Eq. (JT3J). (In the general case, the formula 
()32|) yields the second, linearly independent solution of (|14j)V Conversely, it is easy to 
prove that if any analytic function obeys Eq. (|3"4"j) then Eq. (fT4"j) is also satisfied. 

It is also worth noting the following representation of the phase function through 
the polynomial P n : 

exp(-i^)) = ie^ +1 ^f^. (35) 

It follows from Eqs. (pEjl . (|7jl . (pjl . . 

Eq. (|54~jl leads to the following relations mentioned above among the coefficients of 

Pn- 

e(2o;) _1 ao = — fia n (36) 

e(2^) _1 a fc = (n + 1 - k)a n+1 - k - f^a n -k, k = 1, 2 . . .n, (37) 

where Eq. can be considered as a particular case of Eq. (pT7J). provided one has 
introduced a n+ i = 0. As the above speculation claims, they are equivalent to Eq. (JHJ). 
This implies some further relationships discussed below. 

The elements of the matrix of the linear system (|37 |) . ()36j ) can be represented 
in terms of the Kronecker delta symbols as follows: 

Gfl = e(2ujy 1 5 j k + fiSj n -k - jSj n+i-k, j,k = 0,1, ... ,n. (38) 

Then the product of the (commuting) matrices G^" 1 " 1 -*, G'- -1 -' is easily computable and 
one gets 

(G (+D . G (-D) . l = {X -j(n + l- j))Sj i + (n- j)^ + jfiSj m . (39) 
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The right-hand-side expression here is nothing else but the component representation 
of the matrix <fr (|20p. Thus $ is factorizable, provided the constraint 4uj 2 (X + /i 2 ) = 1 
(see Eqs. 1)1 is taken into account. The same concerns the determinants which are 
factorized as follows: 

A n (A,/i) =detG (+1) -detG^. 

Therefore the spectral equation (j2*U|) is equivalent to the condition 

either det G (+1) = or det G^ 1 ) = 0, (40) 

where the either branch involves the algebraic equation restricting parameters /x, u (and 
depending on n). This property is useful in numerical applications. 

Finally, let us consider polynomial solutions of the two specimens of Eq. ((Hj) with 
common parameter /x but different integer tx's and A's obeying Eq. (J21j) which we denote 
ni, 7X2, ^ 7 ^ 2 ^- An automatic computation establishes the validity of the following 
identity 

[ z - {ni+n ^/ 2 exp(-/i(z + z- 1 )) x (41) 



dz 



HP HP 

J ri2 - 1 rti 2V 1 it>2J6 J- ni J 712 



+ " P P 



n 2 ■ 



where S ni , n2 = z ~( n ^)/ 2 exp(-/x(z + 2" 1 )) x 

[(A« - A< 2 > - i(m - rx 2 )(m + n 2 + 2)) ^ 2 
+i/x(rx 1 -n 2 )^ 1 (l + ^ 2 )] . (42) 

It implies the theorem: 

For /x > i7ie polynomial solutions of Eq. \1$ of different degrees rii,n2 are 
orthogonal on the semi-axis M + witt the weights S„ lin2 J^jj| ) ; i.e. 

'711,512 -fni -^12 ^ ^ ^' (^^) 





Under the same condition the polynomials P n are also normalizable in appropriate norm 
involving the factor exp(— fi(z + z^ 1 )). 

Resuming, the nonlinear first order ODE (Q) arising in the theory of Josephson 
junctions with r.h.s. (J2J) was shown to be equivalent to the linear homogeneous second 
order ODE with polynomial coefficients The latter equation was identified as a 

particular instance of the double confluent Heun equation. It means that the problem 
(QJ), (0) proves completely solvable in terms of the double confluent Heun functions. The 
series of constraints on the problem parameters enumerated by a non-negative integer 
parameter n (see (fTKj) ) was derived whose fulfilment leads to the existence of solutions 
representable in terms of polynomials (P n ) of degree n. The corresponding master 
equation is ((HJ), the constraint equations are (|2"T]) . where A n can be computed by 
means of Eq. (jHUj) through the finite products of 2 x 2 matrices with a single zero 
element, the polynomial coefficients are determined by Eq. (j31j) through the similar 
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matrix products. A curious first order linear "non-classical" two-argument differential 
equation (|34|) which P n has to obey was found. The polynomial solutions of Eq. (JT4"j) 
were shown to constitute the normalizable orthogonal system on the positive semi-axes 
R+. 
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